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Based on the t-t'-t"-J model we have calculated the dynamical spin susceptibilities in the antifer- 
romagnetic (AF) phase for electron-doped cuprates, by use of the slave-boson mean-field theory and 
random phase approximation. Various results for the susceptibilities versus energy and momentum 
have been shown at different dopings. At low energy, except the collective spin- wave mode around 
(tt, n) and 0, we have primarily observed that resonance peaks will appear around (O.Stt, O.Ttt) and 
equivalent points with increasing doping, which are due to the single particle-hole excitations be- 
tween the two AF bands. The peaks are pronounced in the transverse susceptibility but not in the 
longitudinal one. These features are predicted for neutron scattering measurements. 
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Since their discovery the hole- and electron-doped 
high-Tc cuprate superconductors have been noticed to 
exhibit many different properties. This is immedi- 
ately seen from the phase diagram in the tempera- 
ture/doping plane. In hole-doped compounds, e.g., 
La2-a;Sra;Cu04, the antiferromagnetic (AF) order stabi- 
lizes only in a narrow doping region x < 0.02^ whereas 
it persists up to a; = 0.14 in electron-doped ones, e.g., 
Nd2-a:Cea;Cu04i^ The electron-hole asymmetry has also 
been observed in various measurements such as nuclear 
magnetic resonance,'^ inelastic neutron scatteringj-- etc. 
Particularly, recent angle-resolved photoemission spec- 
troscopy (ARPES) experiments^ have revealed the pe- 
culiar Fermi surface (FS) structure in electron-doped 
cuprate Nd2-a:Cea;Cu04. It was found that at low doping 
the FS is a small pocket centered at (7r,0), in contrast 
to the hole-doped casei where it is around (7r/2,7r/2). 
Moreover, upon increased doping another pocket begins 
to form around (7r/2, 7r/2) and eventually at optimal dop- 
ing X — 0.15 the several FS pieces together constitute a 
large curve around (tt, tt). 

Theoretically, the essential role of the next near- 
est neighbor (n.n.) hopping t' has been pointed out 
to understand the electron-hole asymmetryi^ By in- 
clusion of t' and possibly the third n.n. hopping t" 
in two types of models t-J and t-U, numerous nu- 
merical and analytical studies have been carried out 
to explore the different properties of electron-doped 
cupratesy^°i"d^d^d'^d^d^d^ Much attention has been 
paid to the interpretatiorki^iiSiiSiii of the FS evolution 
with doping as observed by ARPESi^ By use of the 
t-t'-t"-U model Kusko et ali^ have derived the mean- 
field (MF) quasiparticle energy bands in the AF state. 
In order to get the consistent results with ARPES, how- 
ever, they need to introduce a doping-dependent effective 
parameter U, see also Ref. Alternatively, we have 

adopted the t-t'-t"-J model to construct the FSfii The 
use of J-type models for the electron-doped cuprates is 
a natural generalization from their extensive application 



to the hole-doped ones, and is largely stimulated by the 
accumulating evidence for the universal d-wave super- 
conducting (SC) gapi^ in both kinds of materials. With- 
out phcnomenological parameters as argued by Kusko 
et al, we have obtained that at low doping only one 
AF band is crossed by the Fermi level around (7r,0), 
and upon increased doping the other one will again be 
crossed around (7r/2, 7r/2). Correspondingly a FS pocket 
forms initially around (tt, 0) and later the new one ap- 
pears around (7r/2,7r/2), in agreement with the ARPES 
data. 

In view of the sample quality, on the other hand, it may 
be questioned whether the multiple pockets around the 
inequivalent points revealed by ARPES are detected from 
the uniform phase of the whole sample or from the dif- 
ferent regions of the inhomogenous sample. Thus other 
experimental measurements complementary to ARPES 
are strongly needed. Naturally, any theoretical predic- 
tion based on the current result will be useful to guide 
the experiments. For this purpose the spin dynamics, 
which can be measured by inelastic neutron scattering, 
is calculated in this paper. If it is really true that the 
new FS emerges in the same uniform phase with the old 
one upon increased doping, the novel particle-hole (p-h) 
excitations between them should lead to characteristic 
features for the spin dynamics which implies the infor- 
mation of the spin/charge excitations. 

At present, the calculations on the dynamical spin sus- 
ceptibilities for electron-doped cupratesiSiiiiS are very 
limited, and mostly done in the SC and normal states. 
Here we will concentrate on the AF phase, and study 
the variation of the susceptibilities with doping which 
has not yet been investigated. The significance of the 
topic is highlighted in the electron-doped materials be- 
cause the AF phase is robust to survive a wide doping 
range and the understanding of its nature becomes cru- 
cial. Thus our motivation is twofold. On one hand, 
we wish to predict some characteristic results for exper- 
imental verification. On the other hand, we perform the 



2 



general formulation of the dynamical spin susceptibilities 
under the background of the AF order, which is rarely 
presented in the literaturei^ The random phase approx- 
imation (RPA) is used to take the spin fluctuation into 
account. Our main result is that with increasing doping 
new resonance peaks will appear at low energy around 
(O.Stt, O.Ttt) and equivalent points, which are pronounced 
in the transverse susceptibility but not in the longitudinal 
one. 

We begin with the t-t'-t"-J model Hamiltonian 
H = -i ^ (cLc.a + h.c.) - i' J2 (clcja+h-c.) 

-m^cla^^a . (1) 

where all the notation is standard. For electron-doping, 
one has t < 0, t' > and t" < oMJML Throughout the 
work |t| is taken as the energy unit. Typical values are 
adopted: t' = 0.3, t" = -0.2 and J = 0.3. 

We treat the Hamiltonian ^ by the slave-boson MF 
theory. Without details we briefly introduce the neces- 
sary formulas. Under assumption of boson condensation 
and definition of MF parameters: the staggered mag- 
netization m — (— l)'(5f) and the uniform bond order 

X = (//ct/jo-) (/■ spinon operator), the Hamiltonian ^ 
is decoupled as follows in momentum spaceii 

Hm = i^kfLfka + ek+qfl+Q^fk+Qa) " 2Jm 

k.a 

Y^^{fLh+Q. + h-c.) + 2NJ{x^ + m") , (2) 

k.a 

where Sk = (2\t\x — Jx)(cosfcj. -I- cosfcj^) — 
4i'a; cos fcx cos fcj, — 2t"x{cos2kx + cos2ky) — ^ {x: 
doping concentration), Q = (7r,7r), and N is the 
total number of lattice sites. means the summa- 

tion over only the magnetic Brillouin zone (MBZ): 
— TT < kx ± ky < n. Above, the local constraint for no 
double occupancy has been treated in average as usually 
done. In addition, we do not consider the potential 
pairing in order to harmonize with the experimental 
observation for the large range of pure AF phase. 



By use of the unitary transformations: fka = 
cos 9kaka + crsin0kPka and fk+Qa = -asmOkaka + 
cos OkPka [o- =T (+) or I (-)], where cos 26*4; = 
(g fc+Q ~ £k)/lk, sin26>fc = -AJm/jk and jk = 

(cfc+Q — £k)^ + (4Jm)2, the Hamiltonian ^ can be 
diagonalized in terms of aka and (3ka, with the two 
energy bands Ck.a = (sfe + £k+Q - 7fc)/2 and ^k^p = 
(efe -HEfe+Q -h7fc)/2. 

For each given doping x and temperature T, the MF 
parameters m and x, as well as the chemical potential fi 
are calculated self-consistently. Then the energy bands 
and the corresponding FS can be plotted in the MBZ. All 
the results have been shown in Ref. TtI for T = 10-3. 
Mainly, with increasing doping the two AF bands become 
close to each other due to the decreasing magnetization. 
Within a very narrow doping range around x = 0.144 
they are both crossed by the Fermi level around (tt, 0) 
(and equivalent points) and (±7r/2, ±7r/2), respectively, 
leading to multiple FS pockets around the inequivalent 
points. In order to reveal their consequence to the observ- 
able physical quantities, we calculate the spin dynamics 
in the following. 

The spin susceptibilities (transverse and longitudinal) 
are defined by 

xt,'^^''\q,q',r) = +l(r.5+W(r)5:(f)(0))(o) , (3) 

where (■■•)(o) means thermal average on the eigen- 
states of i7(M), S+ = E^S+e^"'^^ = Ekfl+.^hl 
, = iS^.V and SI = \T.k.<^fU,Jk.- Cor- 
respondingly the term JJ2{ij) ' rewritten as 

(1/^) Ji<l)iS+SZ^ + S^St^) with J{q) = J(cosq, + 
cos qy). 

We first calculate xt under the quadratic Hamil- 
tonian Hm- By transforming / operator into a and (3 
ones we finally obtain 

5q'..q+QXt^{l^(l + Q^i^n) , (4) 

Xo"'{q,q',ii^n) = Sq, ,qXo" {q,ii^n) , (5) 

where 



X^'{q,iujn) = -^Y^ [cos\ek + ek+q){F^o.+Fpp)+sm\ek + 0k+g){Fap+F0a.)] , (6) 
fc 

1 v-^' 

Xo {q,q + Q,iuJn) = [{sm2ek+q - sm2ek){Faa - Fpf:i) + {sm2ek+q + sm2ek)iFap - F^a)] , (7) 



Xo"{q,i^n) = -r^^'[cos2(0;,-0,,+J(F„„ + f^/3^)+sin2(0fe-0fe+,)(J^„^ + i^^„)] (8) 
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with abbreviations 



F I — 



(77,77' = a,/3) . 



Above LOn are Matsubara frequencies and n{- ■ •) is the 
Fermi function. Since the eigenstates of Hm exhibit the 
AF order that breaks the spin rotational symmetry, the 
simple relation Xo = 2xo^ which holds in the case of 
7M = is no longer valid. Moreover, the nondiagonal cor- 
relation function Xo with q' = q + Q [Eq. 0] arises 
due to the umklapp processes. It is instructive to have 
a look at the physics implied by the diagonal functions 
Xo ^'■^^•'(g,7a;„) from the above Eqs. © and With 
momentum transfer q there are two types of p-h exci- 
tations, either within the single band (a or f3) as de- 
scribed by Faa and Ff^p or between the two bands by 
Fai3 and F^q. We notice that the two, intraband and 



interband p-h excitations, have different weights in their 
contribution to Xo~^^^\q,i^n)- Also, the difference be- 
tween Xo~ i'^n) and Xo^{Qj i^n) exists in their different 
weighing factors for each type of excitations. If we con- 
sider low temperature (T —^ 0) and low energy transfer 
(lu —> 0), the dominant contribution to Xo ^^^\qji^n) 
comes from those excitations within the very vicinity of 
the FS. Thus we expect that new resonance peaks with 
large momentum transfer will arise with increasing dop- 
ing due to the emergence of the new FS pockets around 
(±7r/2,±7r/2). 

We further calculate x'' ^^^^q^q' ^i^n)- The resid- 
ual fluctuation of the J-term beyond MF is consid- 
ered by the RPA. We have obtained the similar equa- 
tions for x^~'"^'^\q,q\i^n) to Eqs. Q) and 101 for 



Xo 



^ ^^^\q,q' ,iu!n), with correspondingly 



X+"(g,7w„) = 

Xo ' (g, iujn) - J{q + Q)[xcl'~ (g + Q, q, »^n)xo ~ (g^ q + Q, »^«) - Xp ~ {q, »^«)X([~ (g + Q, ^^n)] 
[1 + J {q)xo^ {q, iw„)] [i-f J(g -f Q)xo^ (q + Q,iujn)]- J {q)J{q + Q)xt^ {q + Q,q, i^n)xQ^ (q,q + Q, i^n) 



+-( , n ■ ^ Xo {q,q + Q,iuj„) , . 
X [q,q + Q, iuJn) = -. — ^-^ — : jST '(10) 



X^^(g,iw„) = 



denominator of Eq. (jSJ 
XQ^iq,iuJn) 

1 + •^(g)xg^(g,«w„) ' 



(11) 



The formula for x^ (g, i'-^n) becomes complicated due 
to the existence of the nondiagonal xj ■ As a result, 
it may contain a pole, which corresponds to a collec- 
tive spin-wave excitation mode. The imaginary part of 
X^ {q,iujn)\iuj„^uj+ir is shown by Fig.^as a function of 
Lo for various g vectors. Throughout the calculation wc 
have taken T = 10^^\t\, the damping rate F = 0.01 |t|, 
and 500 x 500 k points in the MBZ. Three typical doping 
values are adopted in Fig. ^ small x = 0.04, medium 
X = 0.1 and x = 0.144 where new FS pockets around 
(±7r/2, ±7r/2) appear. Globally it is seen that with in- 
creasing doping the visible susceptibility spans a wider 
energy range, due to the broadening bandwidth and re- 
ducing AF gap. For details, we first look at the region 
q ^ Q — (tt, tt) shown by the several lower panels of 
Figs, n^c) and^d). It is observed that a sharp peak, 
particularly for doping x = 0.04 (note three times reduc- 
tion of the amplitude in this case) , is formed and becomes 
stronger when g is closer to Q. The peak position in each 
panel is around the magnon excitation energy for the cor- 
responding Heisenberg model as shown by the cross on 
the abscissa. This indicates the collective spin-wave ex- 
citation in the presence of carriers. The peak is more 
prominent for x = 0.04 because the collective mode is 
better defined for smaller doping. Similarly, the peak 



from the spin-wave excitation is present at low energy 
when g — ^ 0, as seen from the several upper panels of 
Fig. Ha). 

Then we come to g away from Q and 0. Now the collec- 
tive excitation has a sizable energy. We notice that this 
mode may not always be exhibited by the RPA, for exam- 
ple, in the line: qx+ qy ~ one has x^ = Xo~ because 
of J(g) = 0. However, our interest is the experimentally 
relevant low-energy region where the single p-h excita- 
tions play the essential role when g is away from Q and 0. 
A notable feature observed at low uj (< 0.4 J ^ 40 meV) 
is that the susceptibility becomes finite when x — 0.144 
for some wave vectors, e.g., g = (O.Stt, O.Ttt) as shown in 
Fig. nib) . To fully view the susceptibility in the whole g 
plane, we have fixed the low energy w = O.IJ and plotted 
the density of Imx"' (g, in Fig. [3 There is a huge peak 
around Q in each panel which has been left blank. With 
increasing doping, except an expansion of the blank re- 
gion indicating the broadening of the peak around Q, we 
can observe the considerable intensity around the new 
wave vectors g/27r ~ (0.15,0.37) and equivalent points 
[keep in mind the symmetry of x^ with respect to 
g: qx ^ qy or qx(y) ~qx{y) + Stt]. To more clearly 
see the appearance of the new peaks, we have calculated 
Imx"* (g, oj) vs g along a few lines in the g plane for dif- 
ferent dopings, which are shown by Figs. |nia) andEl^b). 
When X approaches 0.144, the primary feature is that 
new peaks grow rapidly around g/27r ~ (0.14,0.36) and 
(0.36, 0.14) as shown in Fig.OIb). By checking the FS, see 
Fig. 3 in Ref. we conclude that these peaks are in- 



4 




(a) x=0.1 
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FIG. 1: The transverse dynamical susceptibility vs frequency 
for various wave vectors: (a) qy = 0, (b) qx + ly = ti", (c) 
qy — TT, and (d) q^ = qy- The dashed, solid, and red lines in 
each panel represent the results for doping S — 0.04, 0.1, and 
0.144, respectively. Note that all the values for 5 — 0.04 have 
been reduced by three times. The crosses shown in part of 
the panels g and Q denote the single magnon excitation 
energies for the Heisenberg model: Jy^ 4 — (cos q^ + cos qy)^. 
The parameters are: temperature T — 10^'^ and damping 
rate T = 10"^ (in units of 



duced by the interband excitations, e.g, from the pocket 
around (— tt, 0) contributed by the a band to that around 
(— 7r/2,7r/2) contributed by the f3 band. From Fig. [S^a), 
we also see that a hump around q/2n = (0.5, 0) appears 
for X = 0.144, which is due to the intra- /3-band excita- 
tions, e.g., from the FS pocket around {—tt/2, 7r/2) to the 
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FIG. 2: The density plots of Imx^ (9,1^) in the q plane for 
doping X — 0.1 and 0.144, with fixed low energy uj — O.IJ. 
The white region in each panel contains a huge peak around 
(tt, tt) originating from the collective spin-wave excitation, 
which is not plotted in order to highlight the much weaker 
peaks at other wave vectors. With increasing doping, new 
peaks around q/2n ~ (0.15,0.37) and equivalent points are 
clearly seen. 
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FIG. 3: Imx+"('7,'^) [(a) and (b)] and Imx^^((7,t<;) [(c) and 
(d)] at UJ — O.IJ with change of q along the lines: (i) qy = 
[(a) and (c)], (ii) q^ + qy = tv [(b) and (d)]. 



one around {tt/2, 7t/2). Note that all the above new peaks 
are characteristic of the emergence of the new FS pockets 
(or the close proximity of the f3 band to the Fermi level). 
In addition, the peak seen in the region 0.2 < qx/2TT < 0.3 
in Fig. Ola), which slightly moves with doping, comes 
from the intra-a-band excitations between the old pock- 
ets around (7r,0) and (— tt, 0). 

Similarly, we have calculated lmx^^(g,a;) which is 
qualitatively the same as Iuixq'^ {q,uj) but quantitatively 
enhanced in most cases. A few results for the longitu- 
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dinal susceptibility at a fixed energy have been shown 
in Figs. ISJc) andlSJd), for a comparision with the cor- 
responding transverse one. We notice that the primary 
peaks around qjl-K ~ (0.14, 0.36) and (0.36, 0.14) as seen 
in upon increased doping become much weaker in 
X^^, and may even be screened by other peaks nearby. 

Finally, we comment on the new peaks around qjl-n 
(0.14, 0.36) and equivalent points appearing in the trans- 
verse susceptibility. As mentioned above, they arise from 
the interband 'p-h excitations. From the FS plot, we do 
not see a nesting wave vector connecting two FS pockets 
which belong to the two bands, respectively. So the new 
peaks found here are not so strong, but in the same or- 
der as found in some similar calculationsi^i They should 
be observable by neutron scattering which can have mea- 
sured very weak intensity of the spin excitationsnS^ We 
point out that all the plots shown in Fig. O are qualita- 
tively the same for any low energy lo/J< 0.4, but the 
ordinates are enlarged with increasing w. Thus the new 



peaks may be enhanced by a larger uj, making their ob- 
servation easier. 

In conclusion, we have calculated the dynamical spin 
susceptibilities in the AF phase for the electron-doped 
t-t'-t"-J model. Various results for the energy and mo- 
mentum dependences have been given at different dop- 
ings. At low energy, except the collective spin- wave mode 
around (tt, tt) and 0, the primary observation is that 
new resonance peaks will appear around (0.37r, O.Ttt) and 
equivalent points with increasing doping, which are due 
to the interband p-h excitations. These peaks are pro- 
nounced in the transverse susceptibility but not in the 
longitudinal one. We hope that our theoretical results 
will stimulate the experimental measurements. 
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